Spin, charge and single-particle spectral functions of the one-dimensional quarter 

filled Holstein model. 
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We use a recently developed extension of the weak coupling diagrammatic determinantal quan- 
tum Monte Carlo method to investigate the spin, charge and single particle spectral functions of the 
one-dimensional quarter-filled Holstein model with phonon frequency ujo = O.li. As a function of the 
dimensionless electron-phonon coupling we observe a transition from a Luttinger to a Luther-Emery 
liquid with dominant 2kj charge fluctuations. Emphasis is placed on the temperature dependence 
of the single particle spectral function. At high temperatures and in both phases it is well accounted 
for within a self-consistent Born approximation. In the low temperature Luttinger liquid phase we 
observe features which compare favorably with a bosonization approach retaining only forward scat- 
tering. In the Luther-Emery phase, the spectral function at low temperatures shows a quasiparticle 
gap which matches half the spin gap whereas at temperatures above which this quasiparticle gap 
closes, characteristic features of the Luttinger liquid model are apparent. Our results are based on 
lattice simulations on chains up to L=20 for two-particle properties and on CDMFT calculations 
with clusters up to 12 sites for the single-particle spectral function. 

PACS numbers: 71.27.+a, 71.10.-w, 71.10.Fd 
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I. INTRODUCTION 

Including phonon degrees of freedom in model calcu- 
lations of correlated electron systems is challenging but 
necessary for the understanding of many experiments. 
One can mention the quasi one-dimensional organics 
TTF-TCNQ where photoemission experiments are car- 
ried out down to 60 K just above the Peierls transition 
[l| . A detailed modeling of this experimental situation is 
bound to include both electronic correlationsj^, iSj, i^, i^] 
as well as the phonon degrees of freedom p. In two 
dimensions the electron-phonon interactions leads to a 
delicate interplay of superconductivity and charge den- 
sity waves depending on the partial nesting properties 
of the Fermi surface [3, @|- More generally, the ability 
to efficiently include bosonic baths in Quantum Monte 
Carlo (QMC) simulations is a prerequisite for the im- 
plementation of extended dynamical mean-field theories 
(EDMFT) where self-consistency both at the two parti- 
cle (bosonic baths) and single particle levels is required 

The aim of this article is to test on the basis of a non- 
trivial model a recently proposed generalization of the 
weak coupling diagrammatic determinantal QMC algo- 
rithm to include phonon degrees of freedom [ll|, [ij] ■ The 
approach relies on integrating out the phonon degrees of 
freedom at the expense of a retarded interaction and then 
to expand around the non- interacting point. Classes of 
diagrams at a given expansion order can be expressed 
in terms of a determinant, the entries of the matrix be- 
ing the non-interaction Green function. The summation 
over those classes of diagrams is carried out with stochas- 
tic methods. Since the algorithm action based, the CPU 
time scales as {f3L)^ ( /? is the inverse temperature and 
L the number of lattice sites ) and is easily embeded in 
dynamical mean-field self-consistency loops. To obtain 
a full account of the physics, we have carried out lattice 



simulations on lattices up to L = 20 to extract two par- 
ticle quantities and cluster dynamical mean-field theory 
(CDMFT) [IJl calculations on embedded clusters up to 
Lc — 12 to investigate the single particle spectral func- 
tion. As a function of the dimensionless electron-phonon 
coupling and at fixed phonon frequency LOo/t = 0.1, we 
interpret our low temperature results in terms of a transi- 
tion from a Luttinger liquid with gapless spin and charge 
modes to a Luther-Emery liquid with gapful spin and 
gapless charge modes [IJ]. This Peierls phase has dom- 
inant 2kf charge density wave (CDW) correlations. We 
have placed emphasis on the temperature dependence of 
the single-particle spectral function in both phases. At 
high temperatures and in both phases the QMC data 
compares favorably with a self-consistent Born approxi- 
mation 15]. The low temperature properties in the Lut- 
tinger liquid phase compare favorably to a bosonization 
approach retaining only forward scattering [16] whereas 
in the Luther-Emery phase a quasiparticle gap match- 
ing half the spin gap is apparent. The temperature de- 
pendence of the single particle spectral function in the 
Luther-Emery phase is particularly rich; at temperature 
scales where the quasiparticle gap closes, features of the 
Luttinger liquid model are apparent. 



The article is organized as follows. In the next section 
we introduce the model, and briefly review our implemen- 
tation of the CDMFT. We refer the reader to Ref. ^ 
for a detailed description of the QMC method. In Sec- 
tion [Hi] we present our numerical results for two-particle 
and single particle correlation functions across the Peierls 
transition. For completeness sake, two appendices sum- 
marize the self-consistent Born approximation [15] and 
elementary aspects of the the Luttinger model appropri- 
ate for the description of the low-energy excitations of 
the Luttinger liquid phase [la ]. 



II. MODEL AND QUANTUM MONTE CARLO. 



The one-dimensional Holstein model we consider reads: 



fe.cr i i 



with tight binding dispersion relation e{k) = 
—2tcos{ka) — fj,. c\^ creates an electron in Wannier 
state centered on lattice site i and with z-component of 
spin a, cl ^ — -7= J2i e^'^'c] ^ creates an electron in a 

Bloch state with crystal momentum k, hi = X^o-Scr^-'^ 
is the on-site particle number operator and Qi and Pi 
corresponds to the ion displacement and momentum. 

In a recent publication [12l | , we have shown how include 
phonon degrees of freedom in the weak coupling diagram- 
matic determinantal quantum Monte Carlo (DDQMC) 
algorithm [ll|. The key ingredient is to integrate out 
the phonon degrees of freedom at the expense of a re- 
tarded interaction and then to expand around the non- 
interacting limit. We refer the reader to Ref. [I4I for a 
detailed description of the algorithm. 

Since dynamical two particle quantities are notoriously 
hard to compute within cluster methods [l3], we have 
used the DDQMC method to simulate the Holstein model 
on lattices up to L = 20 sites to compute those quanti- 
ties. For the study of the temperature dependence of the 
single particle spectral function, we have found it more 
convenient to adopt the cluster dynamical mean field the- 
ory (CDMFT) on embeded cluster sizes up to Lc = 12. 

CDMFT as opposed to the dynamical cluster approx- 
imation (DCA) is particularly useful to tackle our prob- 
lem. It is a real space method which allows for sponta- 
neous symmetry breaking within a predefined unit cell 
of volume given by the cluster size. To implement the 
method, we decompose the chain into L„, super-cells of 
length Lc- A site, i in the original lattice then corre- 
sponds to a super-cell, R, and an orbital index v running 
from 1 ■ ■ ■ Lc such that: i — R + Ui,. Thereby, the volume 
of the Brillouin zone is reduced by a factor Lc and the 

with 



LcL„ 



quantized wave vectors are given hy K ^ 
n G [—Lu/2, Lu/2[. Within this formulation, the self- 
energy and non-interacting Green function correspond to 
Lc X Lc matrices, i:{K,iuJr„),GQ{K,iuJm)- The CDMFT 
approximation neglects the K dependency of the self- 
energy; 11{K,iuj„i) = S(z(jj„i). In analogy to the DMFT 
approach, one can extract the self-energy by solving on an 
Lc cluster the model at hand subject to a dynamical bath 
QoiiuJm) which has to be determined self-consistently. To 
be more precise: 



Giii^m) = 



1 



Qq {iuJrn) - '^{iu 



r,. 2^ n^l 



Lu ^ Gq (K, iu!„i) - T,{iuj„i) 



(2) 



The last equality corresponds to self-consistency. Hence, 
for a given bath Green function matrix Qo{iuJrn) we use 



the DDQMC method to obtain the corresponding self- 
energy ^(icjm) which in turn, owing to Eq. ^, allows 
us to compute a new bath Green function. This proce- 
dure is repeated till convergence is reached. Within the 
DDQMC the self-consistency is particularly easy to im- 
plement as it is possible to compute the Matsubara Green 
functions directly within the QMC code thus avoiding 
the cumbersome transformation from imaginary time to 
Matsubara frequencies. 

Having determined the self-energy, we compute the lat- 
tice Green functions, g{k,iujm) and A; € [— 7r,7r] with: 



g{k,iuj^ 






-a„) 



M:^— 1 



Go (-R', jtj,„) -i;(itj„j) 



fi^iy 



(3) 



We 



In the above, k = K + m^ with K G 

use CDMFT solely to extract the single particle spectral 
function. The required rotation from the imaginary to 
real time axis is accornplished with a stochastic analytical 
continuation scheme [iSl . Il9| . 



III. NUMERICAL RESULTS 

In this section, we present our numerical results at 
quarter filling, p — 0.5, phonon frequency loq = O.lt, 
which places us in the adiabatic limit, and vary the 
electron-phonon coupling as well as the temperature. We 
first consider spin, charge and pairing correlations as well 
as the optical conductivity and then study in detail the 
temperature dependence of the single particle spectral 
function. Two particle quantities are obtained from sim- 
ulations on an L = 20 site lattice. To at best study single 
particle properties, we have used the CDMFT approxi- 
mation on cluster sizes up to Lc = 12. 

To characterize the strength of the electron-phonon in- 
teraction, we consider the effective mass renormalization 
as obtained from the self-energy diagram shown in Fig. 
M For a flat band of width W, Eq. jMl) yields: 



m* g^ 2 

= 1 + A with A = — — 

TO 2k W 



(4) 



with A the dimensionless electron-phonon coupling. 



A. Spin and charge static and dynamical structure 
factors. 



Equal time charge correlation functions. 



N{q) 






{{hrho) - {hr){ho)) , 



(5) 



are plotted in Fig. [T^. As a function of growing electron- 
phonon coupling, the cusp at 2kp = 7r/2, sig naling a 
power-law decay of the correlation function [20| , evolves 




FIG. 1: Density (a), pairing (b), spin (c) correlation func- 
tions as well as the single particle occupation number (d) as 
a function of electron-phonon coupling. 



towards a clear peak signaling a dominant 2kp charge 
modulation at A = 0.35. Note that at the largest con- 
sidered electron-phonon coupling, a cusp at a higher har- 
monic, 4fci?, is equally apparent. A simple interpretation 
of this charge-density wave stems form the Peierls in- 
stability. For classical phonons the inherent 2k f nesting 
instability of one-dimensional systems renders the metal- 
lic state unstable towards a 2fc/ lattice deformation at 
arbitrarily small electron phonon coupling. In this mean- 
field approach the static lattice deformation triggers the 
opening of a charge gap. It has been argued and shown 
numerically [2l| that this situation cannot be carried over 
to quantum phonons. In this case, quantum fluctuations 
destroy the static lattice deformation and a finite value 



of the electron-phonon coupling is required to destabilize 
the Luttinger liquid. The linear behavior of the charge 
structure factor at long wavelengths (see Fig. [T^) points 
to a metallic state at all considered values of the electron- 
phonon interaction since it amounts to a powerlaw decay 
with modulation q — Q oi the real space charge correlation 
function. 

Since we have not included a Coulomb repulsion in 
our model Hamiltonian, one expects two electrons of op- 
posite spins to share the same lattice deformation and 
thereby bind to form bipolarons. Fig. ^p plots the equal 
time pairing correlation functions in the on-site s-wave 
channel. 



P{r) = (AjAo) with Aj; 



c^ c^ 



(6) 



As apparent, the on-site pairing correlations, P{r — 0), 
grow as the electron-phonon coupling is enhanced from 
A = 0.25 to A = 0.35. This behavior reflects the for- 
mation of bipolarons. On the other hand and in this 
coupling range, the long range pairing correlations are 
suppressed reflecting a tendency towards localization of 
the bipolarons. 

The binding of electrons into spin singlets leads to the 
suppression of the 2fci? spin-spin correlation functions de- 
fined by 



5(g) = 5]e^"-0,.,5. 



z,o; 



(7) 



and plotted in Fig. [TJ:. At A = 0.35 both the q — Q as 
well as the q = 2kp cusps in the spin structure factor 
are smeared out thus lending support to an exponential 
decay of the spin-spin correlation. 

Finally, the single particle occupation number. 



nik) 



cljk,a) 



(8) 



is plotted in Fig. [TJi. As apparent, and on our limited 
lattice size, L = 20, the jump at kp — 7r/4 is dramatically 
suppressed as the electron phonon-interaction grows from 
A = 0.25 to A = 0.35. 

Hence, on the basis of the static quantities, we can 
conclude that a transition between a Luttinger liquid 
metallic phase and a spin gaped CDW state occurs in 
the region 0.25 < A < 0.35. We now provide further sup- 
port for this picture by examining dynamical two-particle 
correlation functions. 

In the Lehmann representation, the dynamical charge 
susceptibility is given by: 



71, m 



-'3^"|(n|njm)|2j(£;„ - ^,„ - u) (9) 



where 



= ^WX]iS*^"'^j ^^'^ ^^6 sum rule N{q) = 



- J dLuN{q,Lu) holds. A similar definition holds for the 
dynamical spin structure factor S{q,uj). 

In the absence of the electron-phonon coupling, both 
spin and charge dynamical structure factors are identical 
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FIG. 2: Intensity plots of the dynamical charge structure fac- 
tor at A = 0.15 (left) and A = 0.35. The x-axis corresponds 
to the momentum q. 



FIG. 3: Intensity plots of the dynamical spin structure factor 
at A = 0.15 (left) and A = 0.35. The x-axis corresponds to 
the momentum q. 



and correspond to the well know particle-hole continuum 
with gapless excitations at q = and q = 2kp. Note 
that at quarter-band filling, 2kp — 7r/2. As apparent 
from the Luttinger hquid model (see Appendix iBJl. the 
phonon mode couples only to the charge degrees of free- 
dom. At weak couplings, A — 0.15, the dynamical charge 
structure factor in Fig.[2]shows precisely this feature; the 
continuum of charge excitations is supplemented by the 
dispersionless phonon mode at uq = O.li. In the spin 
sector (see Fig. |3]) only the continuum of two spinon ex- 
citations is present. 

At larger values of A (A = 0.35) and as a consequence 
of the bipolaron formation spectral weight at low energies 
in the dynamical spin structure factor is suppressed. In 
particular from Fig. [3] we can obtain a rough estimate of 
the spin gap at, A^p ~ 0.2t at A = 0.35. The lattice dis- 
tortion in the Peierls phase is accompanied by a softening 
of the phonon mode. At A = 0.35 (see Fig. [5]) we observe 
a piling up of spectral weight at very low frequencies 
with dominant spectral intensity at q = 2k p. This low 
energy feature corresponds to the slow charge dynamics 
of the bipolaronic 2k f CDW (see Fig. [T]a) 'sdj. The high- 
energy continuum at A = 0.35 in N{q^ lj) is comparable 
to S{q^ uj) at the same coupling. This similarity confirms 
that this structure stems from the particle-hole bubble of 
dressed single particle Green functions. 

We note that phonon dynamics have been studied for 
the spinless Holstein model within a projector based 
renormalization method [22] as well as with exact di- 
agonalization and CPT methods [22]. In analogy to our 
results, the phonon spectral function reveals not only the 
phonon dynamics but also the particle-hole continuum. 

Finally we consider the real part of the optical conduc- 



tivity. 



TT 



n.m 



{l-e-^^)\{n\'j\m)\H{E^^E„,~u) 

(10) 
with i = it Yji,„ (cl,o-Cj+a,cr ~ H.c. j both at A = 0.15 and 
A — 0.35. Our results on an L = 20 lattice are plotted 
in Fig. m As apparent at A = 0.15 a Drude feature 
reflecting polaronic conductivity is visible. In contrast, 
at larger electron-phonon couplings, the formation of the 
bipolaronic CDW leads to a substantial suppression of 
the Drude feature. The suppression of the Drude weight 
reflects the very small charge velocity of the bipolarons. 
This follows from the continuity equation which estab- 
lishes a relation between the optical conductivity and the 
dynamical charge structure factor: 



o'{q,uj) 



1-e 



-I3uj 



)N{q, 



(11) 



At small momentum transfer, and using the sum rule 
^ J dujN{q,uj) = N{q), we can model the dynamical 
charge structure factor by: N{q,u!) = TrN{q)S{vcq — w) 
with Vc the charge velocity. From Fig. [1^ N{q) oc q in the 
long wavelength limit, and the proportionality constant 
is to a good approximation A independent. Inserting this 
approximate form of into Eq. pip gives in the zero tem- 
perature limit: 



lim cr'(g, w) ex VcS{u!). 



(12) 



Hence, the suppression of the Drude weight stems from 
reduction of the charge velocity when passing from the 
Luttinger liquid phase to the bipolaronic CDW phase. 
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FIG. 4: Optical conductivity in the Luttinger liquid phase and 
bipolaronic CDW phases. The calculations were carried out 
with periodic boundary conditions. For this choice of bound- 
ary conditions, the sum rule J dua^uj) = ~Tr{K) where K is 
the kinetic energy, holds only in the thermodynamic limit. In 
the plot, we have imposed this sum rule by normalizing the 
spectra by an overall factor. 



B. Temperature dependence of the single particle 
spectral function 

In this section we study the details of the temperature 
dependence of the single particle spectral function, both 
in the Luttinger and bipolaronic CDW phases. 



1. Atomic limit 

It is instructive to start with the atomic limit, t = 0, 
and in the absence of spin degrees of freedom, 

H^ec^c + gQn+^+'^Q^, (13) 

where exact solutions for the temperature dependence of 
the spectral function are available [2J|. In particular, at 
T = 0, (see Fig. [SJ the single particle spectral function 
is given by: 



AH 



1=0 



[e~A + u„l]) (14) 



with A = g^/2k. 

An electron on the energy level couples to the phonon 
degrees of freedom and can lower its energy at the ex- 
pense of a shift in the ground state expectation value of 
Q. Thereby e — ^ e — A which corresponds to the lowest 
energy pole in A(uj). Since the ground state contains an 
infinite number of phonon excitations, poles at e — A-l-o^o^ 
following a Poisson distribution are apparent in the sin- 
gle particle spectral function. The spectral function is 
centered around (lu) = J dujA{u})w — e and has a width 

d {[to — (lu)] ) = y/uJoA. The relevant energy scale for 
the temperature behavior of the spectral function is the 
phonon frequency, Wq- As apparent in Fig. \E\at temper- 
atures in the vicinity of the phonon frequency a consid- 
erable broadening of the spectral function is apparent. 
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FIG. 5; Spectral function as a function of temperature in the 
atomic limit. 



2. Luttinger Liquid phase, A — 0.25 

Fig. [6] plots the temperature dependence of the sin- 
gle particle spectral function for the Holstein model in 
the Luttinger liquid phase at A = 0.25. We compare our 
results to the self-consistent Born (SCB) approximation 
[15J briefly reviewed in Appendix |^ At high temper- 
atures, T/wq > 1, the overall features of the spectral 
function as obtained from the SCB compare favorably 
with the CDMFT calculations. Both show a broad spec- 
tral function centered around the bare electron energy 
e{k) — e{kp). As in the atomic limit and at an energy 
scale set by the phonon frequency a substantial narrow- 
ing of the spectral function and reordering of spectral 
weight is apparent. 

As the temperature drops well below the phonon fre- 
quency, /3t = 80, the CDMFT spectral function exhibits 
sharp features which are not captured by the SCB ap- 
proximation. For instance at uj/t > and k < kp & sharp 
peak is apparent at w ~ wq in the QMC spectra and is not 
present in the SCB approximation. Of course, the SCB 
approximation has many caveats since i) it does not con- 
tain vertex corrections required in the low-temperature 
Luttinger liquid phase and ii) the phonon propagator is 
not renormalized such that phonon softening and signa- 
tures of the Peierls transitions are not included in the 
approximation. The low temperature CDMFT spectral 
function at A = 0.25 is at best understood within the 
framework of bosonization as sketched in Appendix |B] 
In a first approximation, and deep in the Luttinger liquid 
phase, one can neglect backward scattering ^] thereby 
obtaining the forward scattering model of Eq. (|B10p |16| 
containing spin, phonon and charge modes. The spin 
mode decouples and the charge and phonon mode mix. 
At the expense of a Bogoliubov transformation, the for- 
ward scattering model can be diagonalized to obtain the 
dispersion relations shown in Fig. [7^. Gapless spin and 
polaron modes as well as a gapfull charge mode are appar- 
ent. Since the single electron operator can be expressed 
in terms of the spin and charge operators [26| one expects 
signatures of those modes in the single particle spectral 
function. Fig. [71d plots a closeup of A{k,uo), at our low- 
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FIG. 6: Single particle spectral function at various temper- 
atures. The left hand panels are CDMFT calculations on 
Z/c = 8 clusters. The right hand panels provide a comparison 
with the SCB approximation (See Appendix [Aj The y-axis 
corresponds to the crystal momentum k. 



est temperature. Structures following the coupled gaped 
charge and polaron modes (vertical Hues) are clearly ap- 
parent. According to Fig. [TK the spin mode is next to 
degenerate with the charge modes, and hence difficult to 
detect in our numerical calculations. 
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FIG. 7: Dispersion relation of spin and mixed phonon and 
charge modes as obtained from the the Luttinger liquid for- 
ward scattering Hamiltonian of Eq. [1^. Here we have 
set A = 0.25, LOO ~ O.lf and kp = 7r/4 as appropriate for 
quarter filling. For the purposes of comparison with the 
QMC data, we have taken the liberty of replacing vpk by 
— 2icos(fca) + 2icos(fcFffl) 



3. 



Peierls phase, X — 0.35 



At larger values of the electron-phonon coupling 
backward-scattering becomes relevant and is at the ori- 
gin of the Peierls transition. Fig. [H] tracks the temper- 
ature dependence of the single particle spectral function 
at A = 0.35 which places us in the Peierls phase. At high 
temperatures, 1/(3 > ujq the overall features can again 
be well accounted for within the SCB approximation re- 
viewed in Appendix [X] Upon cooling (see the (31 = 20 
data set in Fig. [5]) a narrow polaronic band crosses the 
Fermi energy, and gaped higher energy excitations show 
precursor features of back-folding. This data set shows 
remarkable similarities with the features observed in the 
Luttinger liquid phase thereby suggesting that aspects 
of the Luttinger liquid spectral functions are apparent 
at finite temperatures above the crossover to the Peierls 
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FIG. 8: Temperature dependence of the single particle spec- 
tral function in the Peierls phase. The results stem from 
CDMFT on an Lc = 12 cluster. The left panels corre- 
spond to logarithmic intensity plots with scale given in Fig. 
|6l The right panels shown the spectral function in a nar- 
row window around the Fermi energy and momentum. Here, 
we have normalized the maximal peak height to unity, and 
the total weight under the spectral function is given by: 



phase. At our lowest temperature, the narrow polaronic 
band develops a gap of the order 2Aqp ~ 0.2i, giving rise 
to rather dispersionless features in the spectral function 
at (jj ~ O.li. We interpret those features in terms of the 
formation of the bipolaronic CDW. Here, removing an 
electron costs the bipolaron binding energy. The fact the 



that the spin gap at A = 0.35 as obtained from Fig. [3] 
matches 2Aqp confirms this interpretation. 



C. Interpretation in terms a transition from a 
Luttinger to a Luther-Emery liquid 

A very natural account of the above presented data 
stems from a transition between Luttinger and Luther- 
Emery liquids. The Luther- Emery liquid description of 
the Peierls phase has been put forward by Voit [ij|. 
Within this framework and away from half filling umk- 
lapp processes leading to a charge gap are absent. Note 
however that at quarter band filling, second order umk- 
lapp processes are allowed and will lead to a charge gap 
provided that the interactions are strong enough such 
that Kp < 1/4 [231 • Here we omit this possibility since it 
does not naturally explain our numerical data on small 
lattices and A < 0.35. Backward scattering on the other 
hand is present and if relevant can lead to the opening of 
a spin gap leaving the charge sector gapless. This corre- 
sponds to the Luther-Emery liquid. 

The Luttinger liquid fix-point is characterized by dom- 
inant forward scattering processes and the asymptotic 
behavior of correlation functions is governed by single 
dimensionless quantity, Kp. Neglecting logarithmic cor- 
rections [28,] the correlation functions read: 



(n(r)n(O)) 



{S{r)Sm = 



K 

— ^ -I- Ai cos{2k fr)r 
(7rr) 

Bi cos{2k fr)r 



-1-K„ 



-1-A' 



(15) 



(A^(r)A(O)) ^Cr-^-^^^" + ■■■ 

Logarithmic corrections do not show up in the first term 
of the charge-charge correlation functions ^2^] and hence 
allow an efficient determination of Kp via: 



-K o = vr iim — z . 

'^ q^o dq 



(16) 



Form our data on an admittedly small lattice, L = 20, 
we obtain from the above equation: 

A' = 1.0341 ±0.0006 at A = 0.15 , , 

(17 
ii:p = 1.0441 ±0.0002 at A = 0.25 ^ ' 

Since Kp > 1 one would conclude that the Luttinger liq- 
uid phase is characterized by dominant superconducting 
correlations. 

The Luther-Emery liquid has correlation functions 
which read: 



An 



(n(r)n(O)) = ^ + Ai cos(2fe/r)r"-^'' 
+A2 cos{Akfr)r^^^'' 
(A^(r)A(O)) ^Cr-^/^" + ■■■ 



(18) 



and an exponential decay of the spin-spin correlations 
|29J . Assuming the validity of the above, we can de- 
duce a rough estimate of the value of Kp in the Luther- 
Emery phase. Since our data at A = 0.35 shows dominant 
2kf charge fluctuations, we conclude that Kp < 1 in the 
Peierls phase. A more precise upper bound for Kp can be 
obtained by comparing the pairing correlation functions 
at A = 0.25 in the Luttinger liquid phase and at A = 0.35. 
At A = 0.25, Kp is slightly larger that unity such that 
the pairing correlations fall of as r~^-^^^ . As apparent 
from Fig. [T|d, the pairing correlations at A = 0.35 in 
the Luther-Emery phase fall off quicker, thus implying 
Kp < 1/2 in the Luther-Emery phase at A = 0.35. This 
upper bound, Kp < 1/2, equally implies a sub-dominant 
Akf charge density decaying more slowly than r~^. The 
observed Akf cusp in the static charge structure factor 
at Akf and A = 0.35 (see Fig. [T^) is consistent with this 
remark. 




FIG. 9: Self-energy diagrams included in the self-consistent 
Born approximation. The solid (wavy) lines corresponds to 
the bare single particle Green function (phonon propagator) 



phase interestingly shows that above the temperature 
scale at which the single gap opens at the Fermi energy, 
features of the Luttinger liquid phase, namely apolaronic 
band crossing the Fermi energy and a gaped charge mode, 
are apparent. This observation should be set in the con- 
text of photoemission experiments carried out on TTF- 
TCNQ organics where measurements are carried out at 
a temperature scale above the Peierls transition and in- 
terpreted in terms of a Luttinger liquid model [l|, y, |3j ■ 



IV. CONCLUSIONS 

In conclusions we have used a generalization of the dia- 
grammatic determinantal QMC algorithm, to investigate 
the physics of the quarter-filled one-dimensional Holstein 
model. We have used the algorithm for lattice simula- 
tions to extract two particle quantities in the context of 
CDMFT to investigate the temperature dependence of 
the single particle spectral function both in the Peierls 
and Luttinger liquid phases. 

Our results are naturally interpreted in terms of a tran- 
sition from Luttinger to Luther-Emery liquids. The Lut- 
tinger liquid phase has a Kp which is marginally greater 
than unity such that pairing correlations are dominant. 
At our considered phonon frequency, loq/I — 0.1, the 
Luther-Emery phase is characterized by Kp < 1/2 and 
thereby by dominant 2kf charge fluctuations. At even 
large values of A than considered in this article, one can 
expect Kp to drop below the 1/4 threshold triggering 
the opening a gap also in the charge sector via second 
order umklapp. Hence at this commensurate filling and 
adiabatic phonon frequency, we can speculate the phase 
diagram as a function of A to not only show transition 
between Luttinger and Luther-Emery liquids but also at 
A > 0.35 a transition from the Luther-Emery phase to 
a fully gaped phase both in the charge and spin sectors. 
One equally expects the character of the Luther-Emery 
phase to very dependent on the phonon frequency. In 
the antiadiabatic limit the Holstein model maps onto the 
attractive Hubbard model where superconducting corre- 
lations are dominant such that Kp > 1. 

Our calculations equally reveal the rich temperature 
dependence of the single particle spectral functions. We 
can access a temperature range covering the domain of 
validity of the self-consistent Born approximation in the 
high temperature limit down to to temperatures where 
the Luttinger liquid or Luther-Emery fix points are rele- 
vant. The temperature dependence in the Luther-Emery 
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APPENDIX A: SELF-CONSISTENT BORN 
APPROXIMATION 

For the Holstein model given by Eq. ([T]), the self- 
energy diagram shown in Fig. [5] can be evaluated to 
give: 



Si(iWm) 



9 
2k''' L 



^E 



nBiuJo) + l-mk)] 



E{k) - UJO + 



UBiLOo) + fHk)] 
-e(fc) +LUO+ ILUm 



(Al) 



Here, f[e{k)] = i3,(l),. is the Fermi function (note that 
we have included the chemical potential in the very def- 



inition of e{k)), nsii-Oo) 



the Bose-Einstein dis- 



tribution and ujQ = \/ -^ the phonon frequency. At zero 

temperature and for real frequencies, the imaginary part 
of the self-energy takes the form: 

2 

ImSi(c^) = -fr^oy J2 {0[£(^)]^(-f(fc) - ^0 + ^) 



2fc L 

k 

+ Q[-e{k)]6{-e{k) + ujo + oj)} . 



(A2) 



The first (second) term in Eq. (|A2[) corresponds to ab- 
sorption (emission) of a phonon. Energy conservation as 
well as phase space limit those processes to energy range 
u; > LOo for absorption and uj < —luq for emission. Hence 
at r = and in a region of width 2luo centered around the 
Fermi energy, the imaginary part of the self-energy van- 
ishes. In this range the single particle Green function has 
poles defining a dispersion relation with effective mass: 



1 



dReE{uj) 



(A3) 



To obtain a good agreement with the high temperature 
Quantum Monte Carlo data we sum up the non-crossing 
self-energy diagrams. This amounts to solving the set of 
self-consistent equations: 



G(fe,za;„i) = 



T.{iujm) 



Here, D{iQrn) 



1 



Gq {k,iuj.„ 



2fc 



I3L ^ 



- S(Jcj,„) 
^ D{in„i)G{k,iujn 



(A4) 



Zi Im)- 



is the bare phonon 



propagator and fi^ a bosonic Matsubara frequency. 
Since at a given iteration we do not have at hand the pole 
structure of G(k, iujm) in the complex frequency plane, it 
is more convenient to solve the above equations numeri- 
cally for real frequencies. To do so, we use the spectral 
representation of the Green function: 



G{k,iu>„ 



duj'- 



A{k,uj') 



ti, LjYi i 



(A5) 



where A{k,uj') = -iG'°'(/s,cj'). With this choice and 
N{uj') = i ^j. A{k,u}') the self-energy reads: 



^r \ 3" I A 'ATI '^ J '^s(^o) + l-/K] 



uja ■ 






(A6) 



At a given iteration step at which N{uj) is known we can 
compute with the above equation the self-energy on the 
real frequency axis {iuJm —* to + id) and thereby recom- 
pute the single particle Green function and corresponding 
N{uj). Typically, for the considered parameter range, ten 
iterations suffice to achieve convergence. 

This approximation has many caveats. Since the 
phonon propagator is not renormalized, phonon soften- 
ing and hence the Peierls transition is absent. At low 
temperatures, in the metallic phase, one equally expects 
the approximation to fail since it does not contain ver- 
tex corrections necessary to produce the Luttinger liquid 
physics. 



APPENDIX B: LUTTINGER-LIQUID 

At low temperatures the self-consistent Born approxi- 
mation does not capture the expected Luttinger behav- 



ior of the one-dimensional Holstein model. Neglecting 
backscattering ~ an approximation which one can justify 
in the Luttinger liquid phase - exact solutions at asymp- 
totically low energy scales are possible. Here, we briefly 
outline the steps. With the bosonic raising and lowering 
operators 



ojoMQi + iPi 



satisfying the bosonic commutation rules 
the Holstein model reads: 



(Bl) 



H= ^<k)^l,a^kM+^oYl 



a\aq + 



Tf Y^ ^k,J'k+q,a (a-l + I 



where the Fourier transform is defined as, 



^k,a 



VL' 



(B2) 



(B3) 



with an equivalent definition for the bosonic phonon op- 
erators flq. 

Linearization around the Fermi points and introducing 
left (Lk.a) and right {Rk,i7 ) fermionic creation operators 
yields the effective low energy form for the kinetic energy 
term, 

k,cr fc.cr 

which in its bosonizcd form reduces to: 
Y'"F\q\bl,a~bq,a with 

q,a 

\l/2 „, 

]f^) J2kRk,aRk+q.a q>0 

, \l/2 ._. . ^^^) 

iffr) J2k^k.a^k+q,a q<0 



(B4) 



bq,a = 



and 



kl 

bq..<r,bl,^, 



^q.q' ^a',cr' 



(B6) 



After linearization the electron-phonon interaction, in 
terms of left and right movers, reads: 



^/2uJoML 



"9+2fe/ 



Yl {^ka^k+q,a [ 
q, fc.cr 

^k.a^k+q,a ( «g-2fcy + ^-q+2kf 



^k.a^k+q,a + ^k.a^k+q.c 



-qr-2fcj 



(B7) 



The first two terms correspond to back-scattering pro- 
cesses which lead to enhanced 2kf charge fluctuations, 
an enhanced effective mass and ultimately to the Peierls 



10 



phase. To obtain a first description of the Luttinger liq- 
uid phase, we omit them thereby obtaining a solvable 
model with only forward scattering processes: 



Hll = ^VF\q\b\^J)q,a + i^o^a\aq 



i-.c 



9 



2LUoMTr \/2 






~q,a ^ "g,cr 



[a.q + a- 



With spin and charge densities defined as, 

Pq = 



bqA - Ki 



1 



K] + bq,l I , 



(B9) 



Hll takes the form: 



(B8) Hll = J^^p'^^^^I^i + J^^^^^^pIpi+^oJ^^I^^ 



9 



2ujr)Mn 



Y.\q\(plq + Pq){al + a^q) (BIO) 



As apparent, the spin mode decouples and the charge 
and phonon modes mix. A Bogoliubov transformation 
diagonalizcs the Hamiltonian and reveals the dispersion 
relation of those modes (see Fig. [T]). 
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